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Abstract 

A part of the two-loop self-energy correction, the so-called P term, is evalu- 
ated numerically for the Is state to all orders in Za. Our calculation, com- 
bined with the previous investigation [S. Mallampalli and J. Sapirstein, Phys. 



CJ 

'^ ' Rev. A 57, 1548 (1998)], yields the total answer for the two-loop self-energy 

>>: 

i-G ' correction in H-like uranium and bismuth. As a result, the major uncertainty 

Ph: 

is eliminated from the theoretical prediction for the Lamb shift in these sys- 
^ ' tems. The total value of the ground-state Lamb shift in H-like uranium is 

2 ■ found to be 463.93(50) eV. 

!>' 

O ■ PACS number(s): 31.30.Jv, 31.10.+Z 

o 

.J:^ ■ INTRODUCTION 

^ . 

P^. The calculation of all two-loop QED diagrams for the Lamb shift of H-like ions is one 



of the most challenging problems in bound-state QED. The experimental accuracy of 1 eV 
aimed at in measurements of the ground-state Lamb shift in H-like uranium [Q] requires a 
calculation of the complete set of QED corrections of the order a^ without any expansion 
in the parameter Za {Z is the nuclear charge number, a is the fine structure constant). In 
high-Z Li-like ions, these diagrams are the source of the major theoretical uncertainty for 
the 2pi/2-2s transition energy [Q] and, therefore, the limiting factor in comparison of theory 
and experiment. Also in the low-Z region, the two-loop Lamb shift is important from the 
experimental point of view 0. What is more, its Za expansion exhibits a rather peculiar 
behavior, with a very slow convergence even in case of hydrogen P] . In order to eliminate the 



uncertainty due to higher-order contributions, it is important to perform a non-perturbative 
(in Za) calculation of two-loop corrections even in the low-Z region. 

The most problematic part of the one-electron a^ contribution is the two-loop self-energy 
correction, represented diagramatically in Fig. ^ The diagram in Fig. |I](a) is usually divided 
into two parts, which are referred to as the irreducible and the reducible contribution. (The 
reducible contribution is defined as a part of this diagram in which intermediate states in 
the spectral decomposition of the middle electron propagator coincide with the initial state.) 
The irreducible contribution (also referred to as the loop-after-loop correction) can be shown 
to be gauge invariant when covariant gauges are used. Its evaluation is not very cumbersome 
and was accomplished in several independent investigations [^j-^. 

The reducible part of the diagram in Fig. |l|(a) should be evaluated together with the 
remaining two diagrams in Fig. ^ This calculation is by far more difficult. The first attempt 
to attack this problem was performed by Mallampalli and Sapirstein [0. In that work, the 
contribution of interest was rearranged into three parts (referred to by the authors as the 
M, P, and F terms), only two of which were actually evaluated. The remaining part (the P 
term) was left out since a new numerical technique had to be developed for its computation. 
In our present investigation, we perform the numerical evaluation of the missing part of the 
two- loop self-energy, the P term. Results of our evaluation, added to those from Ref. 0, 
yield the final answer for the total two-loop self-energy correction for H-like uranium and 
bismuth. This result disagrees with the resent calculation of the total contribution reported 
by Goidenko et al. [^, which is based on the partial- wave renormalization approach. 

The plan of the paper is the following. The basic formulas needed for the evaluation 
of the P term are given in the first section, alongside with the discussion of the treatment 
of ultraviolet and infrared divergences. In the next section we formulate the scheme of the 
numerical evaluation and give some technical details. Numerical results are discussed in the 
last section. In that section we also collect all second-order contributions to the ground-state 
Lamb shift of H-like uranium and to the 2pi/2-2,s transition energy in Li-like uranium. In 
the latter case, the two-loop self-energy contribution is estimated by scaling the Is result. 

In this paper, we use the relativistic units {h = c = m = 1). The roman style (p) is used 
for four vectors, bold face (p) for three vectors and italic style (p) for scalars. Four vectors 



have the form p = (po, p). The scalar product of two four vectors is (p • k) = poko — (p • k). 
We use the notations fi = p^7^, p = p/|p|- 

I. BASIC FORMALISM 

In this paper we are concerned with the evaluation of the correction 

^E = AE^^ + AE^2 ^ 2AE«^ , (1) 

where the contributions AE'^^, AE'^^, and AE'-^^ are represented by Feynman diagrams 
shown in Figs. 0, 0, and ^, respectively. Our consideration of these three sets of Feynman 
diagrams should be considered as an investigation, complementary to Ref. 0, to which we 
refer hereafter as I. In our calculation we use the Feynman gauge and the point nuclear 
model, the same as in I. Our results, combined with those from I, should yield the complete 
answer for the two-loop self-energy correction. 

All the contribution AE^^, AE^'^ , and AE^^ are ultraviolet (UV) divergent. Following 
I, we refer to their UV-finite part as the "P term". We note that subtractions in these 
contributions are chosen in such a way that each of them is free from overlapping divergences. 
The main problem in the evaluation of these diagrams is that they contain bound electron 
propagators as well as UV divergences. While UV divergences are normally separated in 
momentum space, the Dirac-Coulomb Green function is generally treated in the coordinate 
representation. The most direct way for the calculation of the P term consists in developing 
a numerical scheme for the evaluation of the Dirac Coulomb Green function in momentum 
space, which is one of the aims of the present work. 

The nested contributions, AE^^ and AE^"^, possess in addition some infrared (IR) 
divergences, which are associated with the so-called reference-state singularities. These 
divergences are cancelled out when considered together with the corresponding parts of the 
reducible contribution of the diagram in Fig. |l|(a). Following I, we handle the IR divergences 
by introducing a regulator. This makes clear that great care should be taken in order 
to separate all divergences exactly in the same way as in I, in order not to miss a finite 
contribution. 



A. One- loop self- energy 

We start with some basic formulas for the first-order self-energy correction. The formal 
expression for the unrenormalized first-order self-energy matrix element in the Feynman 
gauge is given by 

^E^ren = '^ T dou j d^,d^, i;\{^,) a ^G {s ^ " 0., x^, X2)a^^,(x2) ^^^^^^^^^^ , (2) 

where a^ = (1, ex.) are the Dirac matrices, and G is the Dirac-Coulomb Green function, 

Eq. d^) is written completely in coordinate space. We will need also the corresponding 
expression in the mixed momentum-coordinate representation. This can be obtained by the 
Fourier transformation of Eq. (|^) over one of the radial variables, 

A-(^, q) = , _\ ^ .^ I d^^l{^)a,M^)e-'^- , (5) 

where q = Pi — P2, and 

tPip) = f d^e-'P^'^ix) . (6) 

We note that while the integration over uj in Eq. (^) can be carried out by Cauchy's theorem, 
we prefer to keep it, having in mind future generalizations to the two-loop case. 

The renormalization of the one-loop self-energy is well known. In our work, we employ 
the method based on the expansion of the bound electron propagator in Eq. (^ in terms of 
the interaction with the nuclear Coulomb field . For the detailed description of our renor- 



malization procedure we refer the reader to []I0[. The renormalized self-energy correction is 
represented by the sum of three finite terms, 

where 

^Egl = / ^ V'l(p)7oSg)(e„, p)^„(p) , (8) 



^^^^'^ = I (^ (^ ^iiPihoTUea, Pi; ^a, P2)Vc{^)UP2) , (9) 

where Vc{q) = —ATxaZ/\q\^, and S)j (p) and r^(pi,p2) are the renormahzed free self- 
energy and vertex operators introduced in Appendix ^ The expression for Ai^^^ny i^ given 
by Eq. (Q), where the Green function ^(ea — ^,xi,X2) is replaced with G*^^+^(£:a — co',xi,X2), 

G'(2+)(£,Xi,X2)=G'(£,Xi,X2)-G'(°)(£,Xi,X2)-G'«(£,Xi,X2), (10) 

G'-°^(£:,Xi,X2) is the free Dirac Green function, and 

G(^)(^,xi,X2) = - /rfzG'(°)(£,xi,z)^G(°)(£,z,X2) . (11) 

J z 



B. Basic formulas for two-loop diagrams 

The formal expression for the first diagram in Fig. ^ can be obtained from Eq. (0) by 
the substitution G{ea — u;, xi, X2) -^ GNi{sa — ^, xi, X2) , where 

Gm{e, xi, X2) = / — — ^ G{6, xi, p)7o T.^°\e, p) - 6m G{e, p, xa) , (12) 

G{e, xi, p) = I rfx2 e'^^'^Gie, xi, X2) , (13) 

Gie, p, X2) = I rfxi e-'P^^G{e, x^, X2) , (14) 

and S(°)(p) is the free one-loop self-energy operator defined in Appendix ^. 

The expression for the first diagram in Fig. ^ is obtained from Eq. (|^) by the replacement 

G{ea - uj, xi, X2) -^ GN2{£a " uj , Xi, X2) , wherc 

GN2{,e, xi, X2) = J 7^^?^;^ ^(.^^ ^1' Pi)^c(q)7or°(£:, pi; e, P2)G{e, p2, xa) , (15) 

where r^(pi,p2) is the free one-loop vertex operator defined in Appendix 0, and Vc is the 
Coulomb potential in momentum space. 

We obtain the expression for the first diagram in Fig. ^ by replacing one of the 7 matrices 
in Eq. (^) with the vertex operator F, 



27tJ-oc y (27r)3 (27r)3 V £a - tu - £„(1 - iO) 

xV^n(Pl)7or^(£^a - ^, Pi; £a, P2)^a(P2) , (16) 



with A^"'(u;, q) given by Eq. (||) 

The expressions for the remaining diagrams in Figs. Q-^ are obtained in a similar way, 
by replacing the Dirac-Coulomb Green function G in Eqs. (P^, (]T3|), and ( [T^ ) by (?*■''-' or 

C. Separation of ultraviolet divergences 

In this section we isolate the UV- finite part of AE. Following I, we refer to this contri- 
bution as the P term. Considering the renormalization of the diagrams in Fig. |^, we should 
keep in mind that the inner self-energy loops are always accompanied by the corresponding 
mass counterterms. 

The renormalization of the one-loop self-energy and vertex operators is defined in Ap- 
pendix 0. To handle the UV divergences, we regularize them by working in D = 4 — e 
dimensions. The resulting expressions are 

S(°) (p) -6m = 5(^) {i-m) + T.f (p) , (17) 

r^(pi,P2) = Z(V + r^(pi,P2), (18) 

where B^^^ and U^^ are UV-divergent constants, and S}j (p) and r^(pi,p2) are finite. Ac- 
cording to the Ward identity, B^^^ = —L^^\ 

For the renormalization of the two-loop self-energy diagrams, we refer the reader to a 



(very pedagogical) description given by Fox and Yennie |[TT|]. Applying their arguments to 
the diagrams in Figs. ^|, we have 

AE^' + AE^^ = 5(1) A^ily,^ + AE^' + AE^ , (19) 

2Ai?«i = 2L(i) Ai^il,,^ + 2AE^' , (20) 

where the subscript P means that the corresponding contribution is UV convergent, and 
the subscript D of AE^^^yj^ indicates that this correction should be evaluated in D dimen- 
sions, keeping terms of order e. (These terms yield a finite contribution when multiplied by 
divergent renormalization constants.) 
The resulting expression reads 

AE = L^'^AE^2ny,D + A^p' + A^p' + 2AE^i . (21) 



Here, the correction AEp^ can be obtained from the corresponding expression for AE'^^ 
by the replacement S*^°)(p) -^ S]:j (p), and the corrections AEp"^ and AEp^ - by the cor- 
responding substitution r^(pi,p2) — > r^(pi,p2). Since all the P-terms are UV-convergent, 
in their evaluation we can disregard terms of order e in definition of S]:j (p) and r^(pi, P2). 
We note also that the UV-divergent part of AE separated in Eq. (|2lD , exactly corresponds 
to that in I. 

D. Separation of infrared divergences 

Infrared divergences occur in the corrections AEp^ and AEp"^ due to so-called reference- 
state singularities. They arise when the energy of intermediate states in the spectral de- 
composition of both electron propagators in Eqs. (|T2|), ( P^ coincide with the energy of the 
reference state varCa- As shown, e.g., in I, the divergent terms disappear when considered 
together with the related contributions from the reducible part of the diagram in Fig. |l|(a). 
However, since we are going to evaluate the contributions AEp^ and AEp"^ separately, a 
proper regularization of the IR divergences is needed. In order to preserve the compatibil- 
ity of our results with those of I, we have to employ exactly the same procedure for the 
regularization of the IR divergences. 

Following I, we introduce in the IR-divergent parts of AEp^, AEp"^ a regulator 6 which 
modifies the location of the reference-state pole of the Green function, Sa —>■ £a{^ + ^)- After 
that, we have for the IR-divergent part of AEp^ 

AE-U6) = ^E r du^ l'''^!'^'" ' "^'5 (aa|i^^^e-H---|aa) , (22) 

P,IRK J 27r ^ i-00 {-EaS ~UJ + iOy ^ ' X12 ' ^ ^ ^ 

where a denotes the electron with the energy Sa and the angular-momentum projection /iQ. 
In Appendix ^ we demonstrate that S^ (e^ — u;) as a function of cu can be analytically 
continued to the first quadrant starting from the right half of the real u axis, and to the 
third quadrant from the left half of the real u axis. Therefore, we can perform the Wick 
rotation of the u integration contour in Eq. (^) , 

^E^US) = --T.^ rrf..^^^%-q^(aa|l^^i^e— |aa) . (23) 

' TT -^ Jo {SaS + lUy X12 

Let us investigate the behavior of AEpjp for small values of 6. Writing it in a compact 
form, we have 



/•oo 

Jo 



/M 



{ea6 + iuY 



^f-^fH-m) 



{eaS + iujy 

3fJ[/'(0)] \n6 + terms, regular in 6 . (24) 



Taking into account that 
d 



3? 

dio 



LU=0 



Y, {a\ Sg) {Ea -iuj)\a) {aa\ - — ^i^e"""!^ \aa) 



.(0) 



(a|S^^^(5,)|a) , (25) 



we have 

^E^'miS) = -AEgl In 6 + AE^}^,^ + 0{6) . (26) 

TT 

Here, we have introduced the correction AEpl^^^. that does not depend on the regulator S 
and can be obtained by fitting numerical results for AEp}fj{6). 

An analogous evaluation for the IR-divergent part of AEp"^ yields 



^E^U^) = ^ AE« In 6 + AE^pI,, + 0{6) . (27) 



As can be seen, the IR-divergent parts separated in Eqs. ( p6D and (^) are exactly cancelled 
by the corresponding terms in Eqs. (50) and (55) of I. 

II. NUMERICAL EVALUATION 
A. Green function in the mixed representation 

The main problem of the numerical evaluation of the P terms is that they involve the 
Dirac-Coulomb Green function in momentum space. Until recently, there has been no con- 
venient method for its evaluation. As was pointed out in I, new numerical tools should be 
developed for the calculation of the P terms. 

Here we address two main features which allow us to evaluate the P terms. Firstly, we 
express them in a form which involves the Fourier transform of the Green function over 
only one radial variable [Eqs. (|1^) and ([T^)]. We refer to this as the mixed (momentum- 
coordinate) representation. Secondly, we develop a convenient scheme for the numerical 
evaluation of the Green function in the mixed representation. This scheme was proposed 



and tested in our previous evaluation of the irreducible part of the diagram in Fig. |l|(a) (the 
loop-after- loop contribution) |^. Here we describe the basic idea of this approach. 

We start from the B-splines method for the Dirac equation |T^. For a fixed angular- 



N 



momentum quantum number n, it provides a finite set of radial wave functions \ Lp\ ^(x) \ 
where the superscript i = 1,2 indicates the upper and the lower component of the radial 
wave function, and n numerates the wave functions in the set. The wave functions are found 



as a linear combination of the B-splines |13 



<n(^) = -E«'('^'^'Oi?/(^), (28) 

^ I 

where {Bi{x)}, / = 1, 2, ... is the set of the B-splines defined on the grid {xi}. Since each 

of Bi{x) can be represented as a piecewise polynomial, we can build the corresponding 

piecewise-polynomial representation for the wave functions, 

<n(a;) = - E cU^:, n, I) {x - xi)'' , X e [xi, xi+i] . (29) 

^ k 

Consequently, the radial Dirac-Coulomb Green function in the coordinate space, defined 

as 

G^{e,xi,X2) = 2^^ ^ , (30) 

can be written in an analogous form, 

G'^{e,xi,X2) = ^ A'fc^fc (£,K,/i,/2)(a;i-x/j''i(x2-Xi2)''2 ^ 

xi e [xij^,xi^+i],x2 e [xi^,xi.,+i] . (31) 

Here, the coefficients A*^^^^^ are given by 

At this point, we have built the Dirac-Coulomb Green function in the piecewise-polynomial 
form. This representation is very convenient for the numerical evaluation. After the coeffi- 
cients ^fcjfcj are stored for given values of n and e, the computation of the Green function is 
reduced to the evaluation of a simple polynomial over each of the radial variables. We note 
one additional advantage of this representation of the Green function, as compared to its 
closed analytical form. The Green function in the form (P^) and its derivative are continuous 



functions of the radial arguments, while its analytic representation contains a discontinuous 
function 9{xi — X2) (see, e.g., p^). 

Now we turn to the Green function in the mixed representation. The Fourier transform 
of the radial Green function over the second radial argument can be written as 

roo 

G';l{e,xi,p2) =4:TTs{Lj) J dx2xlJL^{p2X2)G'^{e,xi,X2) , (33) 

where Li^2 = \k± 1/2| — 1/2, s{Li) = 1, 5(^2) = —k/\k\, and Jl{z) denotes the spherical 
Bessel function. Introducing the Fourier-transformed basic polynomials, 

nf (p) = 47rs(Li) [ '^' dxx{x- xifji^px) , (34) 

Jxi 

we write the Green function in the mixed representation, 

G'l{e,Xi,P2) = — $^(Xi -X/J^i 5^ Ag;,^(£, /€,/!, /2)njf^(p2) , Xi G [xi„Xi,+i] . (35) 
^ k\ l2k2 

Certainly, a computation in the mixed representation is essentially more time-consuming 
than that in coordinate space, due to necessity to evaluate the whole set of the integrals 
nJ2^(p2) for each new value of p2- Still, in actual calculations we can perform the numerical 
integration over xi first, and the total amount of computational time can be kept very 
reasonable. 

B. Evaluation of A^^^ 

First, we separate the total contribution AEp^ into two pieces, the IR-divergent part 
AEp]p^ given by Eq. (^21), and the finite remainder AEp^. The expression for AEpl is 
given by 

j^-^ ia f°° J [dp f J J exp{i\uj\xi2) 



x?/;^(^xi)a^^7vi(£:„-u;,xi,p,X2)a''^a(x2) , (36) 



where 



^Ari(£:,Xi,p,X2) = G(£:,Xi,p)7oS2H£:,p)G'(£:, p,X2) 

-G^'\e, xi, p)7oSg^(e, p)G(°)(£, p, X2) 
-G^'\e, xi, p)7oSg)(e, p)G(°)(e, p, X2) 



-G(°)(5,xi,p)7oSg^(e,p)G«(e,P,X2) 

- E M^^lMM 7oSg) (., p) Mp)4(^ . (37) 

The next step is to perform the Wick rotation of the cj-integration contour in Eq. (pG]). 
This is very convenient for the numerical evaluation since, firstly, the Dirac Green function 
as well as the photon propagator are exponentially decaying functions for imaginary values of 
uj. Secondly, by this deformation of the contour we escape most of the problems connected 
with poles of the electron propagators and with the analytic structure of S)j (e). The 
analysis given in Appendix ^ shows that S^ (e^ — cu) as a function of real u allows the 
analytical continuation to the first quadrant of the complex plane starting from the right 
half of the real axis, and to the third quadrant from the left half of the axis. Therefore, 
we can rotate the integration contour on the imaginary axis dividing AEp}. into two pieces, 
AEplj^ corresponding to the integral along the imaginary axis, and the pole term AEp^^i^ 
that arises from the pole of electron propagator with e„ = Sa- (At this moment, we assume 
that a corresponds to the Is-state.) We have 



^'^"^ n Jo J {2-Kf J xx2 



xV^^(xi)a^^jvi(£:„-zu;,xi,p,X2)a''V'a(x2) , (38) 

^^Soie = 2 ^ y 727r)3 y ^^i^^2 — ^l(xi)a/,g^°i^(ea, xi, p, X2)a''^a(x2) , (39) 



where 



^Ar°i''(^a,Xi,p,X2) = ^ [^a(xi)^i(p)7oS2^ (e^, p)G"''='^(£a, P, X2) 



ixeAi -r^\^^ \^y^) 



+G-'^(£„xi,p)7oS^^^(£„p)^^(p)^l(x2)J , (40) 
and G^'^^ is the reduced Dirac-Coulomb Green function, 

G (e,xi,X2)- ^ ,_,^(i_,o)- (41) 
Finally, we collect all contributions to AEp^, 

AE^\6) = -AEgl \n6 + AE^}^,, + AE^^,^ + AE^L + 0{6) . (42) 

We note that instead of dividing AEp^ into three parts, we could have introduced the 

regulator 6 right from the beginning in AEp^, as it was done in I for the "M terms". 



However, this would cause a rapidly varying structure of the integrand in the low-u; region, 
which makes calculations much more time-consuming. {E.g., for the pole term introducing a 
regulator would involve a numerical evaluation of the integral which yields the (5-function in 
the limit 6^0.) In our approach, on the contrary, only the IR-divergent part is evaluated 
with the regulator; the corresponding calculation is relatively simple and allows accurate 
fitting to the form (^61). 

Let us now outline essential features of our numerical evaluation. As can be seen, the 
dependence of the functions Qni and Q^^^^ on the angular parts of xi and X2 is exactly 
the same as that of the Dirac Green function. Therefore, the angular integration causes 
no problems and can be performed by a straightforward generalization of formulas given in 
Ref . |Ty] . The most problematic part of the numerical evaluation of AEp^ is the calculation 
of AEpl^. All numerical integrations were performed by the Gauss-Legendre quadratures. 
The ordering of integrations in our computation coincides with that of Eq. (0), with the 
summation over the angular-momentum quantum number of intermediate states k, moved 
outside. For each value of k and u, we store three sets of complex coefficients ^^^^2 cor- 
responding to the functions G, G^^\ and G^^\ For each value of p, we calculate also a set 
of the Fourier-transformed polynomials 11;^ (p). After this, the integrations over the radial 
variables xi, x^ are performed. This scheme is rather efficient and was used for the evalu- 
ation of AEp^^. The corrections AEpj^^i^^ and AEp^^^^^ were calculated in several different 
ways, which served as a good test for our numerical procedure. 

The actual calculation was performed with the basis set constructed typically with 50- 
60 B-splines of order 6 and employing the point nuclear model. The cavity size of 1 a.u. 



was employed for Z = 80, which was scaled as 7/Z with Z, 7 = wl — (aZy. We use 
an exponential grid with the ffist knot of about 0.001 a.u. for Z = 80. This particular 
grid was chosen since it yields an optimal convergence in the evaluation of the ffist-order 
self-energy correction with respect to the number of knots. The infinite summation over the 
angular-momentum quantum number of intermediate states k, was terminated typically at 
l^tmaxl = 7. The tail of the expansion was estimated by polynomial ffiting in 1/|k|. The 
results of the numerical evaluation of the individual contributions to AEp^ are presented in 
Table | for Z = 83, 90, and 92. The numerical errors, quoted in the table, originate mainly 



from the sensitivity of the result to the number of knots and different grids. 



C. Evaluation of ASP 



The correction AEp"^ can be written in the same way as AEP, 

AEP{6) = -AE« \n6 + AE^l,, + AE^oie + A^^L + 0{6) . (43) 



Here, we again separate the IR-divergent part [Eq. (P?])] and perform the Wick rotation of 
the cu-integration contour, separating the corresponding pole contribution {AEp'^^^^). We 
note that the rotation of the contour is possible because the vertex operator r^{ea—uj, £„— cj) 
as a function of real uj allows an analytic continuation in the region of interest, as shown in 
Appendix 0. The resulting expression reads 

n7V2 "^ f^ J . f dpi dp2 [ j__ j__ exp(-u;a;i2) 



y~^l{yii)a^gN2{ea - iuj, Xi, pi, P2, X2)a^V^a(x2) , (44) 

where 

gN2{£, xi. Pi, P2, X2) = G{e, xi, pi)joT%{e, pi; e, p2)G{e, p2, X2) 

-G^'\e, xi, Pi)7or?j(£, pi; e, P2)G^''\e, p2, X2) 

Assuming that a is the Is state, the pole contribution is given by 

xV'^(xi)a^^^°2''(£:a,Xi,pi,p2,X2)a^V^a(x2) , (46) 



where 



gN°2%ea, Xi, Pi, P2, X2) = Y. [^a(xi)?A|(pi)7or?j(£:„, Pi; Ea, P2)G"''^'^(£^a, P2, X2) 

+G''^{ea,^i,PihoT%{ea,Pi]ea,P2)MP2)ijl{^2)\ . (47) 



The angular integration in these expressions is straightforward, due to the fact that the 
functions Qn2 and Qf^2° have the same dependence on the angular parts of xi and X2 as the 
Dirac Green function G. The numerical evaluation of the correction AEpI^ is much more 



time-consuming than that of AEpj^^. This is because the integration over |p| in Eq. ( pH]) 
(after the angular integration has been carried out) is substituted by the triple integration 
over |pi|, IP2I, and ^ = (pi-p2)- So, the numerical evaluation of Ai^p^^^^ involves one infinite 
partial-wave summation and a seven-fold numerical integration. (One additional integral 
comes from the evaluation of the Green function in the mixed representation.) While it 
is possible to evaluate AEpi^ in a similar way as AEpl^, we have found a more efficient 
method for its computation. To this end, we rewrite Eq. ( ^if ) as follows 



^/v-/ ^ni(Pl)70r?jVn2(p2) , 1 - ttitts _<^xi2 , \ 

^ [Ea -tUJ - e/3) {€a - lUJ ~ E^) Xu J 

where Fp = Tp{Ea — iuJ,pi;Ea — iuJ,p2), ^ni and ■?/'n2 stand for solutions of the Dirac equation 
with the Coulomb potential, m^ and u^ denote solutions of the free Dirac equation, and the 
prime on the sum indicates that the term with e„^ = e„2 = Ea is omitted. In order to 
evaluate Eq. (^8l), we introduce the matrix S, 



S^{^,Pi,P2) = 2^ J — rr-^ : r a^2 e ^^riia), 49 

where V5*(p) (i = 1, 2) stands for the radial components of the corresponding wave function. 
An analogous matrix 5'*^°-' is introduced also for the second part of Eq. (^5]). For each value 
of K and tu, we calculate coefficients of the piecewise-polynomial representation of S*, 5'*^°-'. 
Then, for each values of pi and p2; we store two sets of the Fourier-transformed polynomials, 
n;^(pi) and n;^(p2)- Finally, the integration over ^ is performed. 

The results of the numerical evaluation of the individual contributions to AEp"^ are 
presented in Table || for Z = 83, 90, and 92. 

D. Evaluation of AE'^^ 

The expression for AEp^ can be easily obtained from Eq. (0), after rewriting it in terms 
of the Green function and making the substitutions G -^ G^^^' and F^ -^ Tp, 

AEO^ = -2^a T du f ^^ /dz . ^"P^"^^^) 



oo 



p - ^c^ ^^ "",.,2_^2 



(27r)3 (27r)3 y u^ - <f + iQ 



x?/'^(z)a^G'(^+)(e„-c^,z,pi)7oF^(£:„-u;,pi;£:a,P2)^a(P2) , (50) 



where q = Pi — P2, ^'■^^•'(e, z, pi) is the many-potential Green function [Eq. ([ToD ] in the 
momentum-coordinate representation. The analysis given in Appendix |C| shows that the 
vertex operator T'^{ea — ^,£a) as a function of real u allows the analytical continuation to 
the first quadrant of the complex u plane starting from the right half of the real axis, and 
to the third quadrant from the left half of the axis. Therefore, we can perform the Wick 
rotation of the integration contour, separating the corresponding pole contribution, 

AE^' = AE^^^,,, + AE^l^ . (51) 

The angular integration in Eq. (0) is by far more difficult as compared to the contribu- 
tions considered so far. As the involved expressions are rather lengthy, we do not give their 
detailed consideration here. However, in order to give the reader an idea how the angular 
integration is performed, we note that Eq. (|16D is very similar to the free-vertex contribution 
which is encountered in a calculation of the self-energy screening diagrams (compare with 
Eq. (68) in Ref. |l5l). Basically, the angular integration in Eq. (^) is the same as described 



in detail in Ref. |T5[. The only difference is that in Ref. [|15| the integration was demon- 
strated for two particular states, n = 2pi/2 and 2s. In Eq. (|50|) we need a generalization of 
that procedure for an arbitrary n, which is somewhat tedious but straightforward. 

Let us discuss now the numerical evaluation of AEp^. After the angular integration is 
carried out, a typical contribution to AEp\-^ can be written as follows 



/•oo rco roo rl rco 

t = y^ duj I dpi / dp2 d^ dz 
~^Jo Jo Jo J-1 Jo 



pjplz^ 



up' + g2 
^3i{<iz)^\{AG^^^^^\ea- iu),z,p^)j{p^,p-i,i)ifP^{p-i) , (52) 

where Pi = |pi|, P2 = IP2I, ■C = (Pi ' P2), (f' = Pi + P2 ~ '^PiP2C,, V'* is a radial component 
of the wave function, G|.^+^'^ denotes a radial component of G^'^^^ , Jl is a spherical Bessel 
function, and f{pi,P2,0 originates from the vertex operator. Eq. (^) involves one infinite 
partial-wave summation and a six-fold numerical integration. (The sixth integral is the 
momentum integration in the evaluation of the Green function in the mixed representation. 
Two additional integrations can be also mentioned, one over the Feynman parameter in 
the evaluation of f{pi-,P2-,C)i and another a radial integration in the computation of G^^\ 
This makes the total dimension of the integral to be 8.) Two of these integrations involve 
spherical Bessel functions, which oscillate strongly in the high-momenta region. In order to 



keep the amount of computational time within an acceptable limit, the general scheme of the 
calculation should be chosen carefully. In our approach, we introduce the following change 
of variables |Tg: {pi,P2,^} -^ {x,y,q}, where x = p^+p2, y =p~^-p^, q^ = pI+pI-'^PiP2^, 
p> = nisx{pi,p2), p< = min(pi,p2)- After that, we have 

/•oo /'oo rl roc px px Or i 

/ dpi dp2 d^F{pi,p2,0= dx dy dq- F{pi,p2,0 + FiP2,Pi,0 ■ 

Jo Jo J-1 Jo Jo Jy 2,piP2 L J 

(53) 

In the actual calculation, the outermost loop was the summation over k. The next 
loop is the uj integration. For given values of k and cu, we store a set of coefficients of 
the piecewise-polynomial representation of G']^^^ These coefficients are obtained as the 
difference of the corresponding coefficients for G^, G^k\ ^^^ ^1^^ ^s stated in Eq. (pJ]). The 
next two loops are the integrations over x and y. For each value of pi, we store a set of the 
Fourier-transformed polynomials 11;'^ (pi). The next step is the integration over q and the 
evaluation of f{pi-,P2-,C}i "which involves an integration over the Feynman parameter (see, 
e.g., Ref. [|Tn|)- The innermost loop is the integration over z. Its optimization is the most 
critical part from the point of view of computational time. For small values of g, we use 
Gauss-Legendre quadratures. When q is large, we decompose ji{qx) in a combination of 
sin(gx) and cos(ga;), and use the standard routine for the sin- and cos- Fourier transforms 
based on the generalized Clenshaw-Curtis algorithm. At that stage of the computation, both 
V?^ and GI-^^-*'^ are represented by piecewise polynomials and, therefore, the z integration 
can be performed rather fast. 

A good test for our numerical procedure is to evaluate the many-potential part of the 
first-oder self-energy correction, which can be obtained from Eq. (|50D by the substitution 

The results of the numerical calculation of the individual contributions to AEp^ are 
presented in Table |in|for Z = 83, 90, and 92. 

III. RESULTS AND DISCUSSION 

In this paper we present the numerical evaluation of the correction AE, given by the 
three sets of diagrams which are shown in Figs. ^-^ Putting together Eqs. ([l|), (0), (^2|), 



), and (|51|), we have 



AE = L(^)AEii,., + ^ 






a 

nAf2 , A Z7Af2 , a t:^N2 



ln(5 



+2Ai?gi„i, + 2AEgL + 0(5). (54) 



The UV- finite difference AE — L^^^ AE^ j^ corresponds to what in I is called the P term. 



many 



The IR-divergent contributions, still presented in the P term, are cancelled when considered 
together with Eqs. (50) and (55) of I. When all (5-dependent terms are put together, the 
limit (5 — i> can be taken, and contributions of order 5 and higher vanish. Finite individual 
contributions to AEp^, AEp'^, and AEp^ are listed in Tables |-fT|, respectively. In Table [rV| 
we collect all finite contributions to AE. 

Now we can obtain a finite, gauge-independent (within the covariant gauges) result for 
the sum of the diagrams in Figs. |l](b,c) and the reducible part of the diagram in Fig. |I|(a). 
In order to get this, we should add together the contributions listed in Table |^, the results 
from I for the M terms (Eqs. (50), (52), and (55) of I) and those for the F terms (Table IV 
of I). This yields —0.903(11) eV for the ground state of H-like uranium and —0.575(11) eV 
for bismuth. Adding this to the irreducible part of the diagram in Fig. |l|(a) (—0.971 eV for 
Z = 92 [^,0 and —0.544 eV for Z = 83, this work), we have for the total two-loop self-energy 
correction, given by the diagrams in Fig. |1|, —1.874(11) eV for Z = 92 and —1.119(11) eV 
for Z = 83. 

The numerical results for the sum of the diagrams in Figs. 0(b,c) and the reducible 
part of the diagram in Fig. |l|(a), obtained by combining the present calculation with that 
of I, can be compared with the recent evaluation announced in Ref. [||. The results of 
-0.903(11) eV {Z = 92) and -0.575(11) eV {Z = 83) obtained in this work should be 
compared with 1.28(15) eV and 0.73(9) eV, respectively, reported in Ref. 0. Surprisingly 
enough, the comparison shows that the results disagree even with respect to the overall sign 
of the contribution. Commenting this disagreement, one can mention that the partial- wave 
renormalization procedure, used in Ref. 0, is known to produce certain spurious terms due 
to the noncovariant nature of the regularization (see Ref. [0 and a discussion given in 
Refs. [0,§]). We also note that some assumptions employed in the numerical evaluation of 



Ref. 1^ make it difficult to keep accuracy under proper control in the computation. Still, 
in order to resolve this disagreement, it is desirable to perform an evaluation of the total 
two-loop self-energy correction within the covariant approach from the beginning up to the 
end by the same authors. This will be the aim of our future investigation. 

With this evaluation of the two-loop self-energy correction, we complete the long-lasting 
problem of calculation of all second-order (in a) QED corrections for the hydrogen-like ions 
without an expansion in the parameter Za. The complete set of these corrections is presented 
in Fig. ^. The whole set is conveniently divided into several gauge invariant subsets: SESE 
(a-c), VPVP (d), VPVP (e), VPVP (f), SEVP (g-i), S(VP)E (k). In Table we collect all 
available contributions to the ground-state Lamb shift in 238^91+ rjj^^ nuclear-size correction 
is calculated for the Fermi nuclear model with (r^)^/^ = 5.860(2) fni |T8|. The uncertainty 
of 0.38 eV ascribed to the nuclear-size effect is evaluated as the difference between the 
corrections obtained within the Fermi model and with the homogeneously- charged sphere 
distribution, employing the same rms radius. Some of the a^ QED corrections [VPVP 
(f) and S(VP)E (k)] are evaluated only within the Uehling approximation at present. We 
ascribe the uncertainty of 50% to them. To obtain the binding energy, the Dirac point- 
nucleus eigenvalue of —132279.92(1) eV should be added to the Lamb shift presented in 
Table 0. The error of 0.01 eV of the Dirac binding energy results from the uncertainty of 
the Rydberg constant [^ . As can be seen from the table, the present level of experimental 
precision provides a test of QED effects of first order in a on the level of 5%. 

With the two-loop self-energy calculated for the Is state, we can estimate its value for 
the 2pi/2-2s transition in Li-like ions. As is known, the leading contribution to the self- 
energy arises from small distances, where the self-energy operator is close to a (5-functional 
potential. This gives the well-known 1/n^ scaling for the s-states and zero for the p-states. 
Assuming this scaling, we have a 0.23(20) eV contribution for the one-electron two-loop self- 
energy correction to the 2pi/2-2s transition energy in Li-like uranium. Now we collect all 



second-order QED contributions to this transition energy in uranium, as shown in Table ^ 



In the first line of the table, our previous result of 280.48(11) eV |^ is given, in which all 
available contributions are included, except one-electron second-order QED effects. The 
total value of the transition energy amounts to 280.64(11) (21) eV, which agrees well with 



the experimental result of 280.59(10) eV |3^. In the theoretical prediction, the first quoted 
error arises from the uncertainty due to the finite nuclear size effect and due to higher-order 
electron correlations (see discussion in Ref. 0). The second quoted error corresponds to the 
uncertainty of the second-order one-electron QED effects. 

IV. CONCLUSION 

In this paper we developed a convenient numerical approach to the evaluation of two-loop 
corrections in the mixed momentum-coordinate representation. The elaborated scheme was 
applied to the evaluation of a part of the two-loop self-energy correction which was omitted 
in the previous study |^. It is worth mentioning that our numerical procedure is relatively 
cheap from computational point of view, as compared to |^. While in that work the total 
computational time of about 7000 h was reported for a given value of Z, our evaluation 
requires only 100-150 h on a Pentium-like computer for one value of Z . 

The results of our calculation combined with those of |]^ yield the total value for the 
two-loop self-energy correction of —1.874(11) eV for the ground state of H-like uranium and 
of —1.119(11) eV for bismuth. As this correction has been the last uncalculated second- 
order QED contribution in these systems up to now, our calculation improves significantly 
the accuracy of the theoretical prediction in one-electron ions. The total result for the 
ground-state Lamb shift in H-like uranium amounts to 463.93(50) eV. While the present 
experimental precision of ±13 eV |]1| is not high enough to test the second-order QED 
effects, it is going to be improved by an order of magnitude in the near future []1|. 

The evaluation of the two-loop self-energy correction for the Is state allows us to make 
an estimate of this contribution for the 2pi/2-2s transition energy in Li-like ions. This is of 
particular importance, since in that case the experimental accuracy is much better than for 
H-like ions, which makes Li-like ions very promising for testing second-order QED effects. 
The first estimate of the two-loop self-energy correction allows us to ascribe a well-defined 
uncertainty to the theoretical prediction. For the 2pi/2-2s transition in Li-like uranium, the 
total result amounts to 280.64(24) eV, which should be compared with the experimental 



value of 280.59(10) eV |4 . 
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APPENDIX A: ONE-LOOP SELF-ENERGY AND VERTEX OPERATORS 

The free one-loop self-energy operator in the Feynman gauge is defined by 

^pj ^maj (27r)^k2 + ^0^'^(p-k)2-m2 + 20^ ' ^ ' 

To separate UV divergences, we write the self-energy operator in D = 4 — e dimensions as 

S(o) (p) = 5m + 5(1) {i-m) + T.f (p) . (A2) 

Here, 5m is the mass counterterm, 

5m = ^m(A, + ^) + 0(6), (A3) 

47r 3 

i?(i) is the UV divergent part of the renormalization constant 1 — Z^^, 

5« = -£a, + 0(6), (A4) 

Ag = 2/e — 7s + ln47r, and 7^; is the Euler constant. The contribution S^ (p) is finite; its 
definition agrees with that of I [denoted in I as S^^'°^)(p)]. We note that, evaluating two-loop 
corrections, one should keep terms of order e, since they can yield a finite contribution when 
multiplied by a divergent constant of order 1/e. However, in our present evaluation of the 
P terms we do not need explicit expressions for these contributions. 
The one-loop free vertex operator in the Feynman gauge is given by 

r'^ Pi,P2 = -^^^a / .^^^.2^^n ^- r U2 2 ^ -n ^ r 1^2 2^-n ^ ' ^^ 
J (27r)^ k^ -I- lO (pi — k)^ — m^ + zO (p2 — k)^ — m^ + zO 

We define a finite part of the vertex operator through 

r''(pi,p2) = L(V + ri^(Pi,P2), (A6) 



-1 



where L^^^ is the UV divergent part of the renormahzation constant Zi 

L« = |-A, + 0(e), (A7) 

and the Ward identity B^^' = —L^^' is satisfied. Again, our definition of F^ exactly corre- 
sponds to that of I (A^^) in notations of I). 

The exphcit expressions for the operators S^ , F^ in the hmit e — > can be found in 



Ref. [|TU|. For their exact e-dependent form we refer the reader to F 



APPENDIX B: ANALYTIC PROPERTIES OF THE ONE-LOOP SELF-ENERGY 

OPERATOR 

In this section we consider analytical properties of the one-loop self-energy operator 
^^^\Ea — ^,p) as a function of u. From the definition (|A1|), one can deduce that the 



self-energy operator can be represented by a combination of two basic integrals, 

^ ' ^^ I J (277)^ (k2 + tO) [(p - k)2 -m^ + iO]' ^ ^ 

By introducing the Feynman parametrization of the denominator, we rewrite this as 

^'''^~~rJo ^ J (27r)^ [(k - (1 - x)p)2 - (1 - x){rh? - x^^)f ' 
where m^ = m^ — -iO. Shifting the integration variable k — > k — (1 — a;)p, we obtain 

16vr' f\^ f ^"'k {1,(1 -x)p^} 



U J } = / dx / ——jr ^^^ '-^^ ^ , (B3) 

^ '^^ I Jo i (27r)^[k2-(l-x)(m2_3;p2)]2 ' ^ ^ 



where we have taken into account the identity 

r d^k k 

with A(k^) being a function of k^. The k integration yields 

{J, J,} = ^(4vr)^/^F(l + e/2) £ rfxiMi_|M , (B5) 

where A^ = (1 — x){m'^ — xp^). In the limit e ^ 0, we have 

-iV-^/2 = ^ _ In A^ + 0(e) . (B6) 

Roots of A^ can be easily found. 



N = x{l-x){uj+-uj-iO){uj-+uj-iO) , (B7) 



where uj± = Jm? /x + p^ ± £a . Obviously, ojj^ > u^ = m + Sa and u^ > uj^_ = m — Sa for 
all values of |p| and x G [0, 1]. 

We find that for uj g]— cj^ , u;*J_[, the integrals J, J^ and, therefore, the self-energy operator 
S'-^n^a ^ ^) ^s functions of real uj can be analytically continued both into the upper and 
into the lower half-plane. However, for uo > uo^j^ the self-energy operator allows the analytical 
continuation in the upper half-plane only, and for uj < —uj^_ in the lower half-plane only. We 
can conclude that Tj^^'{ea — oj) is an analytic function of uj in the complex plane with the 
branch cuts [m + Sa — iO, oo — iO) and [—m + £a + ^0, — oo + iO). 

APPENDIX C: ANALYTIC PROPERTIES OF THE VERTEX OPERATOR 

Here we are interested in analytical properties of the vertex operator r'^(pi,p2) as a 
function of uj in two kinematics: a) pi = {Sa — co',pi), p2 = {Sa — ^,P2); and b) pi = 

{Sa-UJ^Vl)^ P2 = (£a,P2)- 

From the definition ([A5| ) , we can deduce that the vertex operator can be written as a 
combination of three basic integrals. 



_ W r ^ {1; k^; k^kj 

''"'''"^ I J (27r)«(k2 + z0)[(pi-k)2-m2 + i0][(p2-k)2-m2 + i0] 



{IJ^lJ^lv} - — ^ / 77TZVn7Dr~~7TTT7I T7^^ '9 ', -mr/^ TT^^ „.9 , -m • (^1) 



Introducing the Feynman parametrization of the denominator and shifting the integration 
variable k ^ k — yq — xp2, we obtain 

^ ' "" ''"^ I Jo ^ Jo ^ J (27r)^ [k2_(^q + a;p2)2_(a;^2_^p2_(^_^)p2)]3 ' 

(C2) 

where q = pi — P2, '"^^ = m'^ — iO, and the identity ( |B4| ) has been taken into account. Shifting 

the integration variable y —>■ xy yields 

// / / \ = ^^f\ [^ dxx f ^"^ {1; ^(^q + P2)A»; a:2(yq + p2)^(yq + p2)^ + k^kj 
^ ' ^' ^"^ % Jo ^ Jo ^"^i (27r)^ [k2-x2(yq + p2)2-x(m2-yp2_(i_^)p2)]3 • 

(C3) 

Now we separate the integral J^j, into two parts, I^y = 1°^, + 1^^,. Here, /^^ corresponds to 
the part of I^i, with k^k,y in the numerator, and /" is the remainder. Integration over k 
yields (L> = 4 - e). 



1 '•I {1; x{yq + ps)^; x^{yq + P2)^(yq + P2)i.} 



^0 



{/././;4 --r(i+e/2)(4.r/- r., r.. ^---^^^-^^-^-;- ^^^^^^^^^^^ , (C4) 



lL = -Til + e/2)iAnr/'^ [' dy [' dx^ , (C5) 



e 2 Jo " Jo N'/^ ' 



where iV = x(t/pi + (1 — y)p2)^ + ^^ — 1/Pi ^ (1 ^ 1/)P2- ^^ ^^e limit e — > 0, we have 

Jj, = ^A, -^-f^dyf^dxx lniV + 0(e) , (C7) 

where A^ = 2/e — 7_b + ln47r. 

Obviously, the denominator A^ is a quadratic polynomial with respect to uo. Let us find 
its roots. In kinematics "a" , we have: 

N = -{l-x){ea-ujf + rh^ + B^, (C8) 

where B"^ = xy{l — y)q^ + (1 — x)ypl + (1 — a;)(l — y)pl, B"^ > 0. We write A^ as 

N = {l-x){uj+-uj-iO){LU_+uj-iO) , (C9) 

where 



m2 + fi2 

o;± = J- ±ea. CIO 

V I — X 

As can easily be seen, uj± > u^. = m ± Sa- 
In kinematics "b" , we have analogously, 

N = y{l- xy){uj+ -u- iO)(cj_ + cu - iO) , (Cll) 

where 



1-^ L2 I ^-"^y 



^± = ^— ^/^^ + -T^-S^K - (1 - ^)^^ + B^)±eA . (C12) 



Again, one can show that uj± > uj^ = m^ea- 

So, we find that V"{ea — uj^Sa — oj) and Vi^Sa — oj, Ea) are analytic functions of oj in the 
complex plane with the branch cuts [m + e^ — iO, 00 — iO) and [—m + e^ + ^0, —00 + iO). 
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TABLES 



TABLE L The individual contributions to AEp, in a.u. 



z 


A£^Kinfr 




A F^l 


Total 


83 


-0.03419 


0.00480 


-0.03951 


-0.06890(20) 


90 


-0.03484 


0.00430 


-0.03780 


-0.06834(20) 


92 


-0.03489 


0.00438 


-0.03737 


-0.06788(20) 



TABLE n. The individual contributions to AEp , in a.u 



z 


A-^P,infr 


A-^P,polc 


^e^:l 


Total 


83 


0.06563 


-0.02127 


0.02178 


0.06614(20) 


90 


0.07200 


-0.02675 


0.02460 


0.06985(20) 


92 


0.07403 


-0.02881 


0.02570 


0.07092(20) 



TABLE in. The individual contributions to AEp, in a.u. 



Z 
83 
90 
92 



A-^P,polc 



^Ep}m 



Total 



-0.03602 
-0.03599 
-0.03605 



0.02496 
0.02417 
0.02392 



-0.01106(15) 
-0.01182(15) 
-0.01213(15) 



TABLE IV. Finite parts of the individual contributions to AE, in a.u. 



z 


AE^^ 


AE^' 


2AE"^ 


Total 


83 


-0.06890 


0.06614 


-0.02212 


-0.02488(40) 


90 


-0.06834 


0.06985 


-0.02364 


-0.02213(40) 


92 


-0.06788 


0.07092 


-0.02426 


-0.02122(40) 



TABLE V. The ground-state Lamb shift in 238u9i+^ [^ gV. 



Finite nuclear size 




198.81(38) 




First-order 


self-energy 


355.05 


m 




vacuum polarization 


-88.60 


m 


Second-order 


SESE (a, irred.) 


-0.97 


i 




SESE (a, red.) (b,c) 


-0.90(1) 


This work + 0] 




VPVP (d) 


-0.22 


[22j2| 




VPVP (e,f) 


-0.75(30) 


\mm 




SEVP (g-i) 


1.12 


m 




S(VP)E (k) 


0.13(6) 


[22J27] 




Total (a-k) 


-1.60(31) 




Nuclear recoil 




0.46 


m 


Nuclear polarization 




-0.20(10) 


[2£-ii 


Lamb shift (theory) 




463.93(50) 




Lamb shift (experiment) 




468.± 13. 


i 



TABLE VL The 2^^ /2-2s transition energy in Li-like ^^^U, in eV. The first quoted error in the 
total theoretical prediction arises from the uncertainty due to the finite nuclear size effect and due 
to higher-order electron correlations. The second quoted error corresponds to the uncertainty of 
the second-order one-electron QED effects. 



Transition energy without second-order one-electron 


QED effects 


280.48(11) 


i 


One-electron second-order 


SESE (a-c) 




0.23(20) 


This work 




VPVP (d) 




0.04 


1111,11 




VPVP (e,f) 




0.10(5) 


H-ii 




SEVP (g-i) 




-0.19 


mm 




S(VP)E (k) 




-0.02(1) 


m 




Total (a-k) 




0.15(21) 




Transition energy (theory) 






280.64(11)(21) 




Transition energy (experiment) 






280.59(10) 


mi 



FIGURES 







a be 

FIG. 1. One-electron self-energy Feynman diagrams of second order in a. 






FIG. 2. Diagramatic representation of the correction AE . For brevity, we do not explicitly 
display the diagrams involving mass conterterms; the inner self-energy loop should be understood 
with the corresponding mass counterterm subtracted. 
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FIG. 3. Diagrammatic representation of the correction AE^"^. 
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FIG. 4. Diagrammatic representation of the correction AE . 
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FIG. 5. One-electron QED corrections of second order in a. 



